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Abstract 

We analyze the viability of the Zee-Babu model as an explanation of observed neutrino masses 
and mixings and the possibility that the model is confirmed or discarded in experiments planned 
for the very close future. The allowed parameter space is studied analytically by using some 
approximations and partial data. Then, a complete scanning of all parameters and constraints is 
performed numerically by using Monte Carlo methods. The cleanest signal of the model will be 
the detection of the doubly charged scalar at the LHC and its correlation with measurements of 
the branching ratio of /i ^ 67 at the MEG experiment. In addition, the model offers interesting 
predictions for t" iJ,^fi^fi~ experiments, lepton-hadron universality tests, the ^13 mixing in 
neutrino oscillations and the {mu)ee parameter of neutrinoless double beta decay. 
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I. INTRODUCTION 

The hints for neutrino masses accumulated in the last 30 years [2|| have been converted 
into a strong evidence in the last 10 years j^, li, [s, ^ [?, [s], lol. ]lo|: the only consistent 
explanation for solar neutrino data, atmospheric neutrino data, reactor and accelerator 
neutrino experiments is based on the hypothesis of massive neutrinos which mix and oscillate. 
Using this hypothesis all neutrino data can be fitted with just two squared mass differences 
and three mixing angles. At present, phases are not needed to explain the data although 
they can be included in the analyses and hopefully will be tested in the future. 

This picture, although certainly a big step forward in our understanding of neutrino 
physics leaves more open questions than the standard model (SM) picture in which neutrinos 
are exactly massless because (1) the SM does not contain righthanded neutrinos, thus, Dirac 
neutrinos are not possible, and (2) the renormalizability of the model, with the minimal Higgs 
content, enforces the exact conservation of the lepton number^ which prevents Majorana 
mass terms. Therefore, to include neutrino masses, we should relax some or several of the 
above ingredients of the SM: 

1. Add righthanded neutrinos. 

2. Add new fields which could allow for lepton number violation while keeping the renor- 
malizability of the model. 

3. Drop the renormalizability of the model. 

The last possibility is very general and allows for neutrino masses without touching the field 
content of the SM. However, its predictivity is very limited and it is only useful when seen 
as a low energy parametrization of a more complete theory containing heavy non-standard 
particles. 

Allowing for righthanded neutrinos seems the simplest and most economical solution: 
if neutrinos are like all the other fermions, it is natural to consider righthanded neutrinos 
coupled to the lefthanded neutrinos and the Higgs scalar, obtaining in this way a Dirac mass 



^ Perturbatively both lepton (L) and baryon (B) numbers are conserved separately. If nonperturbative 
effects arc taken into account only B — L is exactly conserved. Still, L and B violations are tiny at zero 
temperature and density. 
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term as all the other fermions have. One can then easily adjust these Yukawa couplings and 
fit all neutrino masses and mixings without problem. 

However, neutrinos are not like all the other fermions: their masses are extremely small 
rriu < 1 eV. In addition, righthanded neutrinos are completely neutral with respect to the 
SM gauge group, therefore, nothing forbids a righthanded Majorana mass term in the SM. 
In fact, many will argue that it should be written unless a good reason is given to forbid it. If 
this mass is present and if it is very large, one can naturally explain the smallness of neutrino 
masses with the seesaw mechanism and perhaps provide a way to understand the amount of 
matter in the universe through the leptogenesis mechanism. Because of all these virtues, this 
model, SM plus heavy righthanded neutrinos, has almost universally become the standard 
mechanism for neutrino masses. However, the model is not completely satisfactory: the 
heavy Majorana neutrinos contribute, through loop effects, to the Higgs mass with terms 
which are proportional to the square of the righthanded Majorana mass. If this mass is huge, 
the contribution to the Higgs mass will be too large. This is a specific realization of the 
generic hierarchy problem of the SM. However, while within the strict SM one could argue 
that the new heavy particles contributing to the Higgs mass could have very small couplings 
to the Higgs boson, in the seesaw mechanism we have explicit particles with explicit masses 
and couplings giving contributions to the Higgs boson mass and the hierarchy problem 
cannot be swept under the rug anymore. Thus, supersymmetry or any other mechanism 
must be invoked to stabilize the Higgs boson mass. On the other hand, the fact that the 
righthanded neutrino mass is very large makes the effects of righthanded neutrinos negligible 
at low energies, except, of course, those related with neutrino oscillations. The model just 
provides neutrino masses and cannot be tested in present or planned experiments^. 

Adding other fields to the SM gives many possibilities but obviously, if there are no 
righthanded neutrinos, all of them require the non-conservation of the lepton number and 
neutrinos acquire a Majorana mass term. Since supersymmetry provides the best solution 
to the hierarchy problem, one natural choice is to use the fields already present in the 
super symmetric extensions of the SM to generate neutrino masses. Very interesting models 
of this type with spontaneous breakdown of R-parity (see for instance |lll.ll2l|). which implies 



lepton number violation, or explicit R-parity violation (see |13l . Il4l |) have been built. All 
^ Although one can check it indirectly through its effects in leptogenesis. 



3 



these models basically use the sneutrinos to generate neutrino masses. 

Alternatively one could add new ad hoc scalars to the SM in such a way that lepton 
number is not automatically conserved: triplets which develo p a VEV 



singlets which could allow for explicit lepton number violation [17 
of the two mechanisms (see for instance 
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charged 



201 1 or a mixture 



The seesaw mechanism is probably the most natural mechanism for neutrino masses. 
However, the LHC is going to provide results very soon and there are plans to increase 
the precision in /i — 67, /i eee experiments in a couple of orders of magnitude (see for 



instance [2.3]). It may be time to explore alternatives for neutrino mass generation that 
could be confirmed or rejected in planned experiments; among them, the natural candidates 
are models in which neutrino masses are generated through radiative corrections: as the 
masses will be suppressed by loop effects the new particles responsible for them could be 
relatively light and be produced at the LHC, at the International Linear Collider (ILC) and 
may have sizable effects in /i — 67, /x — >■ eee experiments. The simplest model of neutrino 



masses is the Zee-Babu'^ (ZB) model 
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20| which just adds two complex singlet scalar 



fields to the SM (that is, just 4 new degrees of freedom) with neutrino masses generated 



at the two-loop level. Another very interesting model is the Zee model [17|] which adds a 
new scalar doublet and a complex scalar singlet (6 new degrees of freedom), however, the 
simplest version of the model gives a too sharp bimaximal prediction for neutrino mixing 



and has already been excluded 2J, |25|, 1261]. Therefore, we will only consider here the 



Zee-Babu model. In the Zee-Babu model neutrino masses are generated at two loops and 
are proportional to several Yukawa couplings of the new scalars and inversely proportional 
to the square of their masses, therefore the couplings cannot be too small and the scalar 
masses cannot be too large otherwise the generated neutrino masses would be too small. 
This is very interesting because the new scalars may be accessible at the LHC and could 
mediate the processes fi ^ e'y and /i eee with rates measurable in planned experiments. 
In this paper we will sharpen the predictions of the model by using both analytical and 
numerical methods, specially under the assumption that the new scalars are light enough to 
be produced at the LHC. The phenomenology of the Zee-Babu model was recently reviewed 



The model was first proposed in [l^ Q and studied later in 20 1. In the literature, it has been often 
referred to as the Babu model. 
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in [271^ where under certain assumptions analytic limits on several couplings and masses were 
set. The analysis of [^^j is very clever and interesting, however, the calculation of several 
processes in that paper was taken from older papers with some wrong factors of 2 which, 
unfortunately, have also propagated to more recent papers. In addition, ref. makes the 
simplifying assumption that certain couplings are negligible. Thus, we found interesting to 
review the phenomenology of the model by relaxing this assumption. This, unavoidably 
requires a numerical study which will be presented in this paper. We also take into account 
new stronger limits on flavour lepton number violating tau decays from BELLE {29! . Isol and 



B ABAR 3l| . Thus, in section [TTl we sketch the model and review some of its features: the 
neutrino masses, contributions to low energy processes like fi 67, /i ^ 3e, r ^ 3/i and 
favoured values for the parameters of the model. In section [XTTl we study the production and 
decays of the new scalars at the LHC In section [TVl we analyze all the relevant constraints on 
the parameters of the model and obtain its predictions for lepton flavour violating processes 
and for neutrino mass parameters (^13 mixing, Majorana and Dirac phases, neutrinoless 
double beta decay parameter {mu)ee, ■■■)■ Section |V] is devoted to a summary of the results. 
Finally in section |VT] we present our conclusions. 

II. THE ZEE-BABU MODEL 

The Zee-Babu model is the minimal extension of the SM providing neutrino masses and 
mixings compatible with experiment: in addition to the Standard Model field spectrum, it 
only contains one singly charged scalar and one doubly charged scalar. 

In order to fix the notation we will briefly review the Zee-Babu model. We will denote the 
SM particle content as follows: ^ will be the lefthanded lepton doublet, I = 2x2^"^ = ir2Ci'^ is 
just the conjugate lepton doublet used to build Majorana type couplings, e is the righthanded 
lepton and H is the Higgs boson doublet. We also have, of course, the weak gauge bosons W^, 
B^, and the quarks and gluons. As mentioned, the Zee-Babu model contains, in addition, 
two charged singlet scalar fields 

h^, k^^, (1) 
with weak hypercharges ±1 and ±2 respectively. Here we will follow the convention Q = 



See also ref. 



28| 



5 



7^3 + ^ and that h and k destroy negatively charged particles, thus h — h and h'' — 
while k — k and k^ — Then, the Lagrangian can be split into two parts: 

>C = CsM + J^ZB ■ (2) 

The first part, Csu, is the Standard Model Lagrangian: 

CsM = iipi + iep e + {iYe H + h.c.) + • • • (3) 

The dots represent SM gauge boson, Higgs boson and quark kinetic terms, quark Yukawa in- 
teractions and the SM Higgs potential. Generation and SU (2) indices have been suppressed, 
and therefore, F is a completely general 3x3 matrix in generation space. The new terms 
in the Lagrangian are 

CzB = D^h^D^'h + D^k^'D^'k + efeh+ + e^ge k++ + h.c. - Vzb ■ (4) 

Since both h and k are SU (2) singlets the covariant derivative only contains couplings to the 
B gauge boson, which after diagonalization will generate photon and Z-boson couplings with 
the scalars, but no W couphngs. Due to Fermi statistics the Zee type Yukawa couphng, fab, 
is an antisymmetric matrix in flavour space while the Yukawa coupling, gab, of the doubly 

charged scalar k, is a symmetric matrix. The scalar potential Vzb contains all renormalizable 
interactions between the scalars h, k and between them and the standard Higgs doublet: 

Vzb = m'^\h\' + m'^\k\' + Xh\h\' + Xk\k\' + Xh,\h\'\k\' 

+ XuhHH^H + XkH\k\^H^H + {fih^k++ + h.c.) . (5) 

Trilinear terms, like in the Zee model, with two scalar doublets and the h vanish identically 
because the coupling is antisymmetric in SU (2) indices and, therefore, requires two different 
doublets. The last term is particularly interesting because if — > 0, the complete Lagrangian 
has an additional global U{1) symmetry which can be identifled with lepton number L (or 
B — L). In fact, if we assign lepton number 1 to both, the lepton doublet and the righthanded 
lepton singlet, we can also assign lepton number —2 to both, the scalars h and k, in such 
a way that this quantum number is conserved in all the Lagrangian except in the trilinear 
coupling of the scalar potential. Thus, if 7^ 0, lepton number is explicitly broken by the ^- 
coupling in the scalar potential. This is very important because this lepton number violation 
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will be transmitted to the fermionic sector and will finally be responsible for the generation 
of neutrino masses. It is also important to remark that this mechanism for lepton number 
violation requires the simultaneous presence of the four couplings Y, f, g and /i, because 
if any of them vanishes one can always assign quantum numbers in such a way that there 
is a global U{1) symmetry. This means that neutrino masses will require the simultaneous 
presence of the four couplings. 

It is also important to note that the gauge-kinetic part of the Lagrangian is invariant under 
the following U{N) transformations in generation space (for N generations of leptons). 

i^Ud, e^Uee. (6) 

Yukawa couplings, however, break this symmetry. This implies that sets of Yukawa couplings 
related by the following redefinitions are completely equivalent 

(F, /, g) ^ iU}YUe, UjfU,, UjgUe) , (7) 

which in turn means that physical observables should transform correctly under these redef- 
initions. This is important to check the behaviour under flavour transformations of physical 
amplitudes, moreover it can also be used to choose a convenient set of parameters in the 
Yukawa sector and count the number of physical parameters following the methods devel- 
oped in [3^. Thus, using redefinitions of eq. one can choose, without loss of generality, 
Y diagonal with real and positive elements. One could also choose / real and antisymmetric 
and leave g as a. completely general complex symmetric matrix. In addition, one can use 
redefinitions of and fc^"*" to set fi real and positive and to remove one of the phases in g. 
Thus, we finally have^ 12 moduli (3 from Y, 3 from / and 6 from g) and 5 phases (all from 
g) and the real and positive parameter fi (plus, of course, the rest of the parameters in the 
scalar potential). However, we will see later that this convention is not compatible with the 
standard parametrization of neutrino masses and mixings and it will be more convenient 
to use a slightly different convention for Yukawa coupling phases: we will also choose Y 
diagonal with real and positive elements, then we will choose fermion field rephasings to 
remove 3 phases from the elements of gab, leaving the elements of fab complex. Charged 
scalar rephasing can be further used to remove the phase of fi and one of the phases of fab, 

^ The counting can be generalized to n generations of leptons. In that case we will have + n moduli and 
— n — 1 phases. 
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Figure 1: Diagram contributing to the neutrino Majorana mass at two loops. 

for instance we can take /^^ real and positive. Of course the counting of parameters is the 
same as before: we will have 12 moduli (3 from Y, 3 from / and 6 from gat) and 5 phases 
(3 from gab and 2 from fat) and the real and positive parameter /i. 

In any of the discussed conventions, Y is directly related to the masses of charged leptons 
TTLa = YaaV, with V = (H^) = 174 GcV, the VEV of the standard Higgs doublet. Then the 
physical scalar masses are 

= + ^huv"^ , ml = m'^ + XkHV^ . (8) 



A. The neutrino masses 



The first contribution to neutrino masses involving the four relevant couplings appears 



at two loops 



19|, 



201 ] and its Feynman diagram is depicted in fig. [T] 



The calculation of this diagram gives the following mass matrix for the neutrinos (defined 



as an effective term in the Lagrangian C^, = —li'lAiui^L + h.c. 



ujab 



16fifacmcg*Jcd'mdfi 



bd 



(9) 



with 



(27r)4 J (27r)4 (P - ml) - ml) {q^ - mj) {q^ - ml) {k - q)^ - mj 



(10) 



led can be calculated analytically 33j, however, since mc,mrf are the masses of the charged 
leptons, necessarily much lighter than the charged scalars, we can neglect them and obtain 
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a much simpler form 

1 1 n^- 

led ^ (T67r2j2 M^y^^^^ ' ^ = max(mfe,mfe) , (11) 
where /(r) is a function of the ratio of the masses of the scalars r = m\/m\, 

. 1 + 4(log^ r - 1) for r > 1 
7(r) = ^ - ^ ^ , (12) 

1 for r — > 

which is close to one for a wide range of scalar masses. With this approximation the neutrino 
mass matrix can be directly written in terms of the Yukawa coupling matrices, /, g', and Y 

A very important point is that since / is a 3 x 3 antisymmetric matrix, dct / = 0, and 
therefore det A^i, = 0. Thus, at least one of the neutrinos is exactly massless at this order^. 
This is a very important result since it excludes the possibility of degenerate neutrino masses. 

To estimate the value of the largest possible neutrino mass we can take ~ ~ ruh = 
M, then the largest v mass will be 

2 

6.6 X 10-3/'^^ , (14) 

which is the typical seesaw formula, suppressed by some additional couplings and loop 
factors. Because in this model one of the neutrinos is massless, the heaviest neutrino mass 
is fixed by the atmospheric mass difference, thus 777/7/ 0.05 eV and 

f g 150^ > 2 X 10-^ (15) 
mi;. 

since LEP bounds on charged scalar masses are typically M > 100 GeV. This means 
that /'s and gi's cannot be made arbitrarily small and natural values for them can be 
g,f ^ 0.01. Then, for these relatively large couphngs and scalar masses in a range 
M ~ 100 GeV — lOTeV the model will give sizable contributions to low energy processes 
like /X — > eee, — > 67, and scalars that could be produced and detected at the LHC. 



This result does not change if higher orders in charged leptons masses are taken into account in the loop 
integral lab- However, one expects it will change if higher loops are considered. 
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Alternatively, if we assume that the Yukawa couplings are smaller than one, f,g<lwe 
could write 

2 

M ^ 6.6 X IQ-^fg^ < 4 x 10^ TeV , (16) 

which is out of reach of planned experiments. However, constraints on low energy processes 
might require smaller couplings which would lead to much smaller scalar masses. In addition, 
these estimates are very rough; for example the relevant couplings may be related to muon 
physics, and not to tau physics: our estimate on neutrino masses should then be reduced by 
a factor (m^/mT-)^, which of course requires much lighter scalars to match the atmospheric 
neutrino scale. It is therefore very important to carefully establish the parameters and 
validity of the model. 



B. Low energy constraints 

In order to provide neutrino masses compatible with experiment, the Yukawa couplings 
of the charged scalars cannot be too small and their masses cannot be too large. This 
immediately gives rise to a series of flavour lepton number violating processes, as for instance 
fi~ e~7 or fi~ e'^e~e~, with rates which can be, in some cases, at the verge of 
the present experimental limits. This means that we can use these processes to obtain 
information on the parameters of the model and perhaps to conflrm or to exclude the model 
in a close future. In this section we will discuss briefly the relevant processes and collect the 
formulas for our conventions of Yukawa couplings: 



"b c 



: The interesting observable for these processes is the decay width. We 



have (see for instance 



r(C - n^c^l) _ 1 



9abgcd 



Gprnl 



2 



(17) 



r(C-^6"^^) 2(1 + 5,, 

In this expression, the term 6cd takes into account the fact that we may have two 
identical particles in the flnal state. In the case of r decays we have to remem- 
ber that leptonic channels small fraction of the decays BR(£^ — > i'^if. ) = 
R(C ^ K^c^d)^Wa ^l^^) (with BR(/i- e-vv) ^ 100%, BR{t' e-uu 
17.84% and BR(r- fx-uu) ^ 17.36%). 
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Process Experiment (90% CL) Bound (90% CL) 
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Table I: Constraints from tree-level lepton flavour violating decays. 

fi^e^ < — > ^^e^: The k^^ scalar exchange gives also rise to transitions of the type 
H~^e~ — ^ fi^e^ which are well bounded experimentally. The relevant four-fermion effec- 
tive coupling generated by exchange of the scalar k~^~^ is (here we use the conventions 



for the effective Hamiltonian and the limits of 



35 



36|) 



^ _ y/2gee9l^ /ION 
^MM o ■ [^°) 



ml 



We collect the relevant constraints of this type in table HI 

—>■ ^bi^v: These processes receive additional contributions from the exchange of 
the singly charged scalar which affect the Fermi muon decay constant but do not 



2l| 



modify the spectrum 

where a sum over undetected neutrinos has been performed. The second term is clearly 
subdominant if nih ^ 200 GeV, however we have included it in the numerical analysis 
and have checked that we can neglect it in analytical estimates. In this model the 
charged scalar only contributes to lepton decays but does not contribute to hadronic 
decays, therefore the effective Gp extracted from hadronic decays and G^ are different. 
However, in the framework of the SM, the equality of Gp and has been tested with 
very good accuracy once all radiative corrections have been correctly included. In the 
SM, both iVudp + \ Vus\^ + iKibP = 1 and Gpsu = G^sm are satisfied. Thus, assuming 
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Gf3SM = Gf^sM one can test the unitarity of the CKM matrix, or, conversely, assuming 
the unitarity of the CKM matrix one can test the universahty of couphngs in hadronic 
and leptonic decays. In the model we are considering this is not true anymore; the 
CKM matrix is still unitary but, as explained, G/3 = G^sm = G^sm 7^ G^. Since the 
extraction of the experimental values of |VJj^^P assumes the SM, we will have 

yexp ^^y,, 

where Vij are the truly unitary CKM matrix elements in the model. Therefore 



(20) 



ICT + IKTI' + IKTT 



n2 



G2 



1 



(21) 



and, since \VlY? + + \Kh^? = 0-9992 ± 0.0011 [37] is very close to 1, we will 

obtain a strong limit on |/e^p/m^. 

On the other hand, the charged scalar contribution will also modify the Fermi coupling 
extracted from r decays in the different leptonic channels. After subtracting the 
different factors from phase space and radiative corrections this is usually expressed 
in terms of ratios of effective "gauge couplings" g^^^ for the different leptons which in 
the SM are all equal (see for instance [38|). Thus, comparing tau decays to muons and 
tau decays to electrons we have (since in the SM Ga^b oc gagt) 



Similarly 



„exp 

exp 
9e 



exp 



G 



Gr- 



1 + 



V2 



G 



1 + 



V2 



( I //ir I I fer \ ) 



exp 



exp 



G 



G 



/i— >e 



1 



Gprnl 
Gpml 



( I fer I I fefj, \ ) 
(l/^trl ~ Ife/il ) 



(22) 

(23) 
(24) 



Universality constraints are summarized in table [TTl , where measured values are translated 
into 90% CL limits. 



— >■ £^ 7 : In the case of transition amplitudes a ^ b the interesting observable is the 
decay rate. We find (for calculations including singly and doubly charged scalars see 



for instance 



21 



39 



m 



C7) 



487r 



(/V) 



ah 



Gprni 



16 



(9^9) 



ab 



Gprnl 



(25) 
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SM Test Experiment Bound (90%CL) 

lept./hadr. univ. E<?=d,s,fe I^T^P = 0.9992 ± 0.0011 |/e;,p < 0.015 {mh/TeVf 
lije universality g""^"^ / gT^ = 1.0001 ± 0.0020 || - < 0.05 {mh/TeVf 

t/^jl universality g'f^lgt^ = 1.0004 ± 0.0022 - < 0.06 (m/,/TeV)2 

r/e universality g'f^lgf^ = 1.0004 ± 0.0023 \\f^,r\^ - |/e^p| < 0.06 (mft/TeV)^ 

Table II: Constraints from universality of charged currents. 
The factor 16 in front of the doubly charged contribution does not usually appear 



in the literature 27| and deserves a comment: the Feynman rule for the kea^h vertex 
contains a factor 2 when a^h because there are two identical terms in the Lagrangian, 
but also the vertex with a = h contains a factor 2 because there are two identical Wick 
contractions in this kind of vertices. This factor of 2 for identical particles was missed 
in which led the authors do define new coupling constants with different factors 
of 2 for diagonal and non-diagonal terms^. It is also important to remark that the 
singly and doubly charged scalar contributions do not interfere because they couple to 
fermions with different chirality. Again we have to remember that BR(£~ ^(Tt) — 

[I — e conversion in nuclei: The new scalars of the model do not couple to quarks and, 
therefore, do not generate a four-fermion operator that could contribute at tree level 
to fi — e conversion. However, radiative corrections, in particular those related with the 
fi — e — '-f vertex, will contribute to the process. It is also clear that those corrections are 
tightly related to the /i — 67 decay discussed above but are not identical because the 



photon in yU — e conversion is not on the mass shell. In fact in ref. 4l|] if was shown that 
in models with doubly charged scalars there is a logarithmic enhancement, log(g^/mfc), 
of the /i — e conversion amplitude with respect to the /i — 67 amplitude. At present the 
best limits come from fi — e conversion on Ti, a{fi~Ti —>■ e~Ti) / a{e~Ti —>■ capture) < 



4.3 X 10~^^ 37(1 which, when translated into limits on the couplings, are slightly worse 
than present yU — > 67 constraints, but one has to keep in mind that if yU — 67 is relevant 
and if /i — e conversion limits are improved in the future it will also be relevant. 



^ One can also see tliat tliose results cannot be right because pliysical amplitudes should transform correctly 
under the flavour redefinitions of couplings in eq. 
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Experiment 



Bound (90%CL) 



6ae = (12 ± 10) X 10-12 r (|/,^|2 + \f,,f) + 4 (l^.eP + IffeMp + \9er\^) < 5-5 X 103 (mfc/TeV)^ 



5a, 



(21 ± 10) X 10-10 r (|/e^|2 + |/^,|2) + 4 (|5e;.P + + l^^rP) < 7.9 (mfc/TeV)^ 



S/?(/i ^ 67) < 1.2 X 10-11 r^|/,Wr + 16|5e*e5eM + 5e*M5MM + ^er^Mrl' < 3.4 X IQ-^ (mfc/TeV) 
BR{t ^ 67) < 1.1 X 10-7 rVe*M VP + 16|<7:e5er + 9*e^9^.r + <7:r5rrP < 1.7 (mfc/TeV)^ 

Bi?(r ^ ^7) < 4.5 X 10-8 r2|/e*u/erP + lQ\9*e.9er + 5* 5/.r + 9*9rr\^ < 0.7 (mfe/TeV)^ 



Table III: Constraints from lepton number violating photon interactions. 

a = {g — 2)/2: For diagonal transitions the model gives additional contributions to 
the anomalous magnetic moments of the leptons which are very well measured in the 
case of the electron and the muon. We find that the additional contribution to the aa 
of lepton ia, in this model, is 



5aa = - < 



SM 



247r2 



{ff)aa , Ag^g) 



+ 4- 



mr 



mr 



(26) 



It is important to remark that we always find a negative contribution (this is in agree- 



ment with old and well tested calculations 
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431]). In the case of the muon a^, recent 



analyses of experimental data and theoretical calculations in the Standard Model sug- 
gest that the exp erimental measurement is slightly larger than the SM prediction (for 



a review see 



44(1). Several authors have tried to explain this la to 3cr effect in differ- 



ent extensions of the SM. In particular, in [45] the charged scalars of the Zee model 
were used to increase the of the SM. We find this is not possible and instead we 
will use the g — 2 oi the muon (and also of the electron) to constrain the parameters 
of the model. The relevant constraints coming from flavour lepton number changing 
photon interactions are summarized in table IIIII . Experimental limits on branching 
ratios are already provided in the literature jsT] at 90% CL. Results for g — 2 are 



usually given as measurements on 5aa. We will use 



46] 5ae = (12 ± 10) x lO'i^ and 



Sa,. 



(21 ± 10) X 10 1°. Notice that in both cases the central value is positive, while 



the model gives a negative contribution. To place 90% CL bounds in this situation 
we use the Feldman and Cousins prescription 4^] which, for the values above, gives 
\Sae\ < 6.1 X 10-12 and \Sa^\ < 3.8 x 10~i°. 

Since lepton number is not conserved, another interesting low energy process that could arise 
in the model is neutrinoless double beta decay (Oz/2/5). In this model, the singly and doubly 
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charged scalars do not couple to hadrons, this means that the (Oi/2/3) rate is dominated 
by the Majorana neutrino mass exchange and it is proportional to the |(A^i/)eeP matrix 
element, therefore it will be addressed in section HVl when we discuss in detail the neutrino 
mass matrix constraints. 



C. Pert ur bat ivity constraints 

Beside requiring that the model produce acceptable (a) neutrino masses and mixings (b) 
low energy predictions, we also have to address theoretical questions related to the validity 
of the predictions and the consistency of the model. Indeed to be able to perform any 
calculation in this model we have to assume that perturbation theory can be used. This 
imposes strong constraints on the relevant couplings of the model. The Yukawa couplings 
of the new scalars receive loop corrections like 

^3 3 

Sf , 5q ^ — — . (27) 

"J (47r)2' ^ (47r)2 ^ > 

If the corrections are going to be much smaller than the couplings, the couplings must 
satisfy /, (? -C 4-7?. Similarly, the trilinear coupling among charged scalars proportional to 
the parameter induces loop corrections to the charged scalar masses like 

^^l^^<-^2- (28) 
Requiring that the corrections are much smaller than the masses implies /i <C 4nmh,4:7rmk. 
Since it is difficult to establish the exact values of the couplings for which the perturbativity 
of the theory breaks down, we will encode this type of constraints in the single parameter k 
and will require 

|/ab|<«, \gab\<n, fj.< Kmm{mh,mk) . (29) 

For the purposes of illustration we will take k = 1 or k = 5, the value k = 5 being rather 
conservative (for instance, the strong coupling constant, Qs, is considered to become non- 
perturbative at scales of about 1 GeV, at those scales as{l GeV) ~ 0.5 and gs{l GeV) ~ 2.5). 

The parameter /i is important in the generation of neutrino masses, therefore the con- 
straint /i < K,mm{mh,mk) is important. On the other hand, as seen in tables [H [Tll IIIIl if 
the scalar masses are relatively light (around 1 TeV or less), low energy processes already 
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provide interesting limits on the charged scalar Yukawa couplings. However, if the charged 
scalars are heavier, the experimental limits on the new Yukawa couplings are so mild that 
may allow Yukawa couplings large enough to compromise the perturbative validity of the 
theory. Then, the perturbativity constraints we just discussed will become relevant. 



III. THE MODEL AT THE LHC 



Extra scalar degrees of freedom arise in many scenarios extending the weak interactions 
beyond the SM. In our case, the scalar sector is enlarged by the addition of two charged 
scalars: h and k, which could be produced at the LHC if their masses are low enough. In 
particular, as we will see below, the LHC will be very well suited for searching the doubly 



charged scalar, k. Studies 



directed in the past [48, 



49 



'or searching doubly-charged scalars at future colliders have been 



50 



51 



52| . In general this scalar is taken to be one component 



of a weak triplet. Such triplets are well motivated on theoretical grounds, specially when 
considering left-right symmetric models. Our work differs essentially in the gauge charges of 
the scalars. Both, h and k, are charged weak singlets that do not acquire a VEV. This makes 
the phenomenology different. More model independent studies have been also considered in 



the literature 



53 



5^. 



Concerning experimental bounds, LEP searched for these scalars. Their pair production 



55 



56 



57| . Single production 



(e+e" 7*Z* — > k^^k ) implies the bound > 100 GeV 

via e~^e~ — ^ kee as well as contributions to Bhabha scattering have been also studied by 



LEP 



57|, 



but in these cases the bounds depend on the (unknown) values of the Yukawa 



couplings. Tevatron has also been used to set bounds on this kind of scalars [59|, |60!, l61 |. 
Depending on the details of the model (couplings, decay channels,...) the mass is again 
found to be roughly above 100 GeV. 



A. Collider phenomenology 



1. Production 

The extra scalars can be pair produced via a Drell-Yan process, fig. [21 Although this 
production mechanism presents the drawback of having a potentially high threshold due to 
the creation of two scalars, it has the important advantage of being proportional to their 
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k+- 



Figure 2: Pair production of k 

gauge charges as well as depending only on one unknown parameter: the mass of the scalar. 
The partonic cross section at LO reads 



a = 



2n2^3 r2Ql _ 2{g, + gn)Q, s - Ml {gl + gl 



(30) 

where s is the energy squared in the center of mass frame (CM) of the quarks, Q stands for 
electric charges, gi and gR are given for the quarks hy gL = T3 — s^Qg and gR = —s'^Qq 



and (3 is the velocity of the produced scalars in this frame /5 = yl — 4m^/s. 

Equation (130|) shows that pair production is four times more efficient for k than for h due 
to their charges (assuming equal masses), which translates into a better discovery potential 
for k. The k pair production cross section, akk, at NLO for the LHC and Tevatron is displayed 
in fig. [31 To compute it, we have used CompHEP [62] with CTEQ6.1L libraries [63] to find 
the LO cross section and afterwards we have included a K- factor of 1.25 for the LHC and 
1.3 for Tevatron to take into account NLO corrections, see 64]. 



Single production might be also interesting when double production is not possible. Single 
production can proceed with a k accompanied by two singly charged scalars, fig. HI or by 
two charged leptons replacing the scalar /I's. If the k is accompanied by two charged leptons 
the amplitudes are proportional to the Yukawa couplings, whose exact values we ignore and 
might be small. 

It is important to note that the cross section will be dominated by the virtual particles 
in the propagators if they could be on-shell. In the case of k being produced with two h, the 
single production will be dominated by the first diagram if s > 2mk, because in this case 
k* can be created on-shell. One might argue that the energy in the center of mass frame of 
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Figure 3: Pair production cross section for k. We have used CompHEP (CTEQ6.1L) to obtain the 
LO and appUed a K-factor of 1.25 for the LHC and 1.3 for Tevatron. 




Figure 4: Single production diagrams. 

the colliding quarks is not fixed, instead it is a fraction of the total energy in the center of 
mass frame of the colliding protons, s. However, the cross section involves an integration 
over the possible values of s. If s is large enough to create two /c's, the integration will be 
dominated by the real production of two /c's, thus reducing the single production to pair 
production. Specifically, a{k'^'^h'~h~) ~ akkBr{k hh). The same reasoning is valid in the 
case of single production with leptons. We have performed calculations using CompHEP to 
check this point. Therefore, single production is only important when the available energy, 
s, is not sufficient to create a pair of k. 
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A possible third production mechanism is via the couphngs with the Higgs doublet, H. 
There is little information concerning these couplings because their contribution to low 
energy phenomenology is expected to be negligible in front of the Yukawas fab and gab- 
This is so because the former enter at two loops and the later at tree level. In any case, 
the amplitudes of these processes are expected to be small because the Higgs couplings to 
quarks are proportional to their masses. 

In summary, we find that, from the point of view of production, the best suited channel 
for discovery studies is pair production, being four times more efficient for k than for h. 



2. Decay 

We assume for the moment that the scalars are not long lived, i.e. they decay before 
reaching the detector. Different decay channels present different experimental sensitivities 
depending on the final products. In particular, detectors are much less sensitive to those 
channels containing neutrinos and/or taus in the final states, since neutrinos (including 
those coming from the decay of the taus) will escape undetected. This makes necessary to 
compute the branching ratios. 

The k scalar can always decay to two leptons of the same sign, since > 100 GeV. The 
width reads 

It is worth to stress that k is the only particle in this model that can decay to two hke-sign 
leptons, which will be crucial to detect it. 

If ruk > ^ruh, then k can also decay to a pair of h 



r{k -^hh)^ — 

Stt 













rUk \ 


_mk_ 





1-^. (32) 



On the other hand, the decay channels of h reduce to those with one lepton and one 
neutrino in the final state: 

V{h ^ ^ayb) = ^m^. (33) 

Since these channels involve always one neutrino it is clear that detecting h will be much 
more complicated than detecting k even if their production rates were similar. 
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Finally we will check that the doubly charged scalar k^~^ cannot be long-lived in this 
model^ once low energy phenomenology and neutrino data are taken into account. Indeed 
the k decay width can be written as 

r = ^(W= + ---)- (34) 

where the gf^^ term takes into account the decay into muons and the dots represent all other 
possible couplings. Then, the long-lived condition translates into Ig^fil < 10^^ which cannot 
be fulfilled when the limits obtained in tabs. [rv¥Vl are used. 



3. Detection at colliders 



As we have seen the discovery potential of the LHC for k is more promising than for h. 
On one side because the production cross section of the former is enhanced with respect to 
the later and, on the other side, because the experimental sensitivity to the decay channels 
of h is smaller. Thus, in the following we will focus on k. 

After a pair of fc's is created in the collider, they can decay into a number of final states. 
The most interesting for us contains four like-sign leptons. From now on we will refer to this 
channel as Alep. The rest of the possible final states always contain h or r. These channels 
are quite difficult to deal with experimentally because r and h will decay to neutrinos. In 
contrast, detection of electrons and muons is quite efficient. In addition, the decay of a 
to a pair of like-sign leptons (e^ or /i^) with high invariant masses constitutes a clear and 
distinct signature. This channel has a negligible background coming from SM processes, 
making it very appropriate for k discovery studies. 

In order to model the efficiencies and acceptances of the detectors at the LHC we use the 
criterion that 10 events of k pair production with subsequent decay to 4lep lead to discovery 
of k. We expect that at least two of such events will be properly detected/identified providing 
us with four pairs of like-sign leptons which invariant mass will give us the first estimate of 



ruk- This criterion is taken from 50[ where the authors perform a study of the discovery 



reach at Tevatron for doubly- charged bosons decaying to like-sign leptons in a similar model 



and a similar rule can be extracted from 



52| where the authors focus on the ATLAS detector 



By "long-lived" we understand that the scalar can travel a distance of the order ct > 3m 61 1 
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Figure 5: Number of events at the LHC in the 4:lep channel for a luminosity C = 300 fb and 
^/s = 14 TeV assuming that all produced k pairs decay in this particular channel. 

at the LHC. A more detailed study of the forthcoming detectors acceptances and efficiencies 
at the LHC is desirable. 

To estimate the maximum reach in terms of at the LHC we take the most optimistic 
scenario in which all the k pairs decay to 4lep. The number of events in this channel is 
shown in fig. [5] for the optimistic luminosity^ 300 fb~^ and CM energy 14 TeV. From this 
plot one concludes that the LHC will be able to probe masses up to 1 TeV approximately. 

In general, the signal in the 4lep channel will be smaller than the one shown in fig. [5] due 
to the presence of the other decay channels. This signal draining will be controlled by the 
branching ratio BR4if.p, which can be expressed in terms of the couplings as 



\9ee 




\9ee 


^ + \9,f.\' + 2\ge,\' + \gHh 


^ + 2\ger 


' + 2\g^r 


2 + 





(35) 



where we have defined the effective coupling of the doubly charged scalar to singly charged 



^ The LHC luminosity is expected to be about 100 fb ^/year. 
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scalars, ghh, as 



( 



Am 



2\ 1/4 



dhh — 




(36) 



IV. ANALYSIS OF THE PARAMETERS OF THE MODEL 

As discussed in section [TTl the Yukawa couplings of the model can be written in terms 
of 12 moduli and 5 phases. Other parameters relevant for the model are the masses of 
the charged scalars, rrih and m^, and the coupling /i. Among the 12 moduli, 3 correspond 
to the 3 charged lepton masses, which are known. Thus, we have 17 additional relevant 
parameters with respect to the plain SM with massless neutrinos (9 moduli and 5 phases 
from the new Yukawa couplings, the 3 scalar parameters, rrih and and the coupling /i). 
The remaining parameters in the scalar potential are of no interest for our purposes. The 
neutrino mass matrix is rather well known. In our case it contains 2 neutrino masses, 3 
real mixing angles and 2 phases (1 CKM-type phase and 1 Majorana phase). Thus, there 
will still remain 4 moduli and 3 additional phases in the Yukawa couplings (plus rrih, 
and /i) . On most of these Yukawa couplings we have some information from section III Bl 
as long as the masses of the charged scalars are not much heavier than 1 TeV. Notice that 
this is the interesting range for scalar masses if they are going to be produced at the LHC 
In addition, there are also indirect arguments that suggest that the scalar masses should 
be relatively light (below 1 TeV) if one likes to avoid strong hierarchy problems^'', since 
the charged scalar masses will contribute, at one loop, to the mass of the SM Higgs boson. 
However, the couplings of the SM Higgs boson to the new scalars are unknown and could 
be small. Thus, although the natural range of the masses of the new scalars is about few 
TeV or less, they can also be larger. Then, in what follows, we will allow the masses of 
the charged scalar to vary between the LEP lower bound ~ 100 GeV and infinity. We will 
immediately see, however, that present information already constrains the charged scalar 
masses to be below ~ 10^ TeV. The couplings gab, fab and must in addition satisfy the 
perturbativity constraints discussed in section III CI 

From the previous discussion it is clear that even though we have 17 additional parameters 
we also have a lot of information on them both from neutrino oscillations and from low energy 

Alternatively one could enlarge the model by supersymmetrizing it. 
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processes. However, the correlations among the different observables due to their dependence 
on the same set of parameters can be difficuh to disentangle. Under those circumstances an 
adequate approach to the analysis of the Zee-Babu model should involve both an analytic 
understanding and a systematic numerical inspection of the parameter space to clarify the 
ranges allowed by the available experimental results for as many parameters and observables 
as possible. 

We first exploit the information we have on neutrino oscillation experiments: the knowl- 
edge of two squared mass differences and the mixing angles. 

With antisymmetric /ab's, at two loops, the mass matrix determinant is equal to zero, 
and thus one eigenvalue is zero. Two mass differences are then sufficient to fix the masses 
both in the normal hierarchy and in the inverted hierarchy cases (a degenerate spectrum 
cannot arise in this model). Except for the Dirac phase 6 and the Majorana phase 0, we 
can almost reconstruct experimentally the neutrino mass matrix by using the information 
we have on the mixing angles and the masses. Without loss of generality we can write the 
neutrino Majorana mass matrix as 



(37) 



with U the standard PMNS matrix 



U 



1 

C23 S23 
-S23 C23 J 



( 



ci3 si3e-^^^ ( 



1 







V 



Cl2 Si2 
-Sl2 Ci2 

1 



(38) 



y -5136*'' Ci3 j 

while Dy is the diagonal matrix of masses (including the only Majorana phase). Notice that 
writing the mass matrix in this form already implies some phase convention. 

Since one of the v masses of the model is zero we only have two possibilities^^: 

Normal hierarchy (NH) 



D 



NH 





m2e"t' 
m3 



, m3 > m2 



A5 
Aa 



ml 



(39) 



Here we follow the conventions and results of ref. 



6a | adapted to our case. 
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Inverted hierarchy (IH) 



mi 

















A5 = m\ — mf 
, mi Ri m2 , . (40) 

= ml 



y u u u y 



With 



As = (7.9 ±0.3) X 10-5 eV^ 
Aa = (2.5 ±0.25) X 10-3 eV^ 



(41) 



s?2 = sin^ ^12 = 0.30 ±0.03 

s?3 = W < 0.02 , 90%CL . (42) 
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= sin^ 023 = 0.50 ± 0.08 
Thus, apart from the poorly known S13 mixing (we just know it is small) and the phases, 
6 and 0, the mass matrix can be partially reconstructed in terms of the 2 known mass 
differences and the 2 known mixing angles for each of the two cases. In particular, we can 
immediately extract the matrix element responsible for (Oz/2/9) decays: 
Normal hierarchy 

NH\ / A aNH\ /a~„2 „2 JS , /a „2 -i25 



(m^^),e = (-Mr )ee = Ascfs^f^e^'^ ± ^/^Asi.e-''' . (43) 

In this case^^, given the previous values, it is clear that {Ai^^)ee ~ 0.003 eV and therefore 
difficult to see in (Oz/2/3) decay experiments. 
Inverted hierarchy 

(m^^),, = {M'fU = VAa ± Ascl,sl2e"^ + V^^^ • (44) 

In this case, unless a cancellation occurs between the two terms for e**^ = —1, {Ail^)ee is 
naturally of order 0.05 eV and, therefore, observable in planned (Oz/2/?) decay experiments. 

Equation ( fT3|) gives the mass matrix Ai,^ in terms of the parameters of the model - the 
Yukawa couplings, the scalar masses and the trilinear coupling -, we can thus try to fix some 
parameters by matching the A4u, obtained from the neutrino oscillation parameters, to the 
calculated one. Since the mass matrix is symmetric, in principle this gives 6 equations. 
However, one of them is trivially satisfied because, by construction, both matrices already 



Notice the dependence on the Dirac phase S. This is a consequence of our convention for Majorana phases. 
One could redefine phases and make this quantity independent on 5, but this wiU not affect predictions 
or constraints on observables. 
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satisfy det{M.i,) = 0. To choose the remaining 5 equations we will use that the eigenvector 
corresponding to the eigenvalue is very simple; as det / = 0, there is an eigenvector a of / 
with zero eigenvalue /-a = 0, a = (/^^, —fer, fe^i)- Obviously, a will also be an eigenvector 
of Ail, with zero eigenvalue when expressed in terms of masses and mixings and, therefore, 
UD^U^a = or, since U is unitary, D^U^ai = 0. This gives us three equations, one of 
which is satisfied trivially because one of the diagonal values of D,y is zero. The other two 
equations will allow us to express the ratios of couplings fij just in terms of mixing angles 
and phases^^. Thus, in the NH case we have {D^^)ii = and 

(f/^a) =0 ^ ^ = tanel2^^ + tan0l3sin023e-^^ 

JfiT COS 6^13 

(U^si) = ^ ^ = tan - tan ^13 cos ^236"*^ • (45) 

COS6'i3 

These equations immediately tell us that the standard PMNS convention of phases is not 
compatible with all fab being real. However, we can take a phase convention in which /^,- is 
real and positive and leave fer and fe^ complex with phases fixed by eq. ( H5l) . 

With values like sfg ~ 0.3, S23 ~ 0.5 and sf^ < 0.02, the first term on the right hand side 
of eqs. (l45l) dominates and we get fer — f^irl'^ — fe^i- With this relation we can go back 
to the low energy bounds in table [TTl and table IIIII and find that the strongest constraints 
on the fij couplings come from /i — > 67 (which strongly bounds |/e,-/^^|) and tell us that 
\fei\ < 0.05(m;,/TeV) and \f^r\ ^ 0.1(m;,/TeV). 

The equations corresponding to the inverted hierarchy case, {Dl^)33 = 0, are 

In this case, it is clear that fer/feiM = -tan 6^23 ~ -1 and |/e^| > 5|/^^|, \fer\ > 5|/^r|. Now 
we can use these relations in the low energy bounds in table [III and table [Inl and find that the 
strongest constraints on the fij couplings come from lepton-hadron universality (see table 
nil), 51/^.1 <|/e.|<0.1(mjTeV). 

We still have 3 additional equations we can use to fix the parameters of the model. In 



Therefore, the decay branching ratios of the scalar h to the different leptons are fixed by the mixing angles. 
This can probably be exploited 28| at the ILC . 
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this case we have no good argument to choose them^^ and following 27|] we will take the 
three elements 17122, ^23 and 77133 in the equalities niij = {J^u)ij = C fiaT^aQab'^bf jb^'^^^ the 
rriij defined by eq. ([ST]) and C = 43^2^/2 1- Thus, if ujah = magl^^nib we have 



"^23 — C{~fnT^tJ.T — ftirfer^eT + ftirfe/i^eii + ferft 



^33 = C(/^r^M/^ + "^fl^r ferret, + fer^ee) ■ (47) 

Because of the hierarchy among the charged lepton masses, it is natural to assume that 
those Uab containing the electron mass, c^ee, ^^e^i, ^^er, are much smaller than u^^, u^^, Urr, 
in that case we can neglect them (we will check later the goodness of this approximation 
within the numerical analysis), and we have 

In the normal hierarchy case this gives {sij = sin^j^, Cij = cos 9ij) 

Cfl^UJrr ^ '^3Ci3S23 + m2e*'^(Ci2C23 - e"^Si2Si3S23)^ , 

Cflr^l,T - -m3c\^C23S23 + m2e*'^(ci2S23 + e"'c23Si2Si3) (C12C23 " e"'si2Si3S23) , 
Cflr^^,^, ^ mscl^^cls + 17126"^ {C12S23 + e'^C23Sl2Sl3f. (49) 



With 7713 ~ 0.05 eV and 7772 ^ 0.009 eV, 

0.05 eV 



(50) 



setting a definite hierarchy among the Qab couplings: 

9tt ■■ 9f,T ■■ fi-MM ~ ml/ml : m^/nir : 1 . (51) 
In the inverted hierarchy case, eqs. give 

Cflr^TT - mi{c23Sl2 + e'^Ci2Si3S23f +m2e^'^{Ci2C23- e'^Si2Si3S23f , 

Cf'^r^^.r - "^l(si2S23 - e*''ci2C23Si3)(c23Si2 + e*'^Ci2Si3S23) 

+ ?772e*'^(Ci2S23 + e*''c23Si2Si3)(Ci2C23 - e"'si2Si3S23) , (52) 

(f^r^^f, ~ mi(si2S23 - e*'^Ci2C23Si3)^ + ?772e*'^(Ci2S23 + 6**^023512513)^ 



Except that they cannot be in the same column or the same row of Al^, because in that case the equations 
are related by Mi^a — 0. 
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where mi ~ m2 — 0.05 eV, also yielding for e*"^ = 1 



0.05 eV 



and the hierarchy of couplings in eq. (!5T|) . However, in the IH case there is a strong cancel- 
lation for Majorana phases close to vr, and one can obtain a smaller value for u;^^, thus we 
can only write 

In both cases one expects (7^^ to be the largest coupling among the three considered. Of 
course, g^e, geij. and g^r can also be large and are only constrained by low energy processes and 
perturbativity constraints. One should notice, however, that in the inverted hierarchy case, 
the approximation made in going from eqs. P7|) to eqs. fHHj) may be a priori less justifiable 
than in the normal hierarchy case when ^13 ^0, as the eigenvector corresponding to the zero 
eigenvalue, (/^^, -fer, fe^i) is proportional to (e*''tan6'i3, sin 6^23, cos6'23), i-e. f^r oc tan6'i3, 
and since the terms retained in eq. (1471) are proportional to f^r-, it is not obvious that the 
terms proportional to uj^i can be neglected. 

Assuming then that ^ \g ^^li'm ^ / rriT-) and \gTT\ ~ Ifi'/x^tK'^^t/'^r)^, we can go back to 
tables m and IIIII to find the relevant constraints on the couplings. The best constraint comes 
from — fi^fi~fi~, which tells us that 1(7^^1 ^ 0.4(mA;/TeV), \gfj,r\ ^ 0.024(mjfc/TeV), 
\grr\ ^ 0.0015(mfc/TeV). 

We can use all this information to set analytical bounds on the relevant parameters of 
the model in the line discussed at the beginning of section III A[ 

A. Analytical constraints 

1. NH case 

First, just from the neutrino mass formula, we have 

,, „ ,9 omaxfrn^., m/,) maxfrn^., m»,) , . 

\9,,\\Ur\' > 10-' 

Now we can show that due to the logarithmic growth of / for ^ rrih and the fact that 
/ < 1 for nik < mh, 
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Thus 



i9..ii/.j>io-^g^ (57) 



Now we use the perturbativity bound on /x, < K,mm{mh, mk] 



^M\U>io-'^'^^^!^^^, (58) 

lev min(mfc, rrih) 

which can be rewritten as (use that mhrrik = m.a.x{mk, mh) min(mfc, mh)) 

/max(mfc,m,,)y 3 |2"^fc /^^x 
J < 10 n\g,,\\f^,\ ^ (59) 

We can use that < max(mfc,m/i) and the perturbative constraints \g^^\ < k, If^irl < 
to find an upper hmit on the masses of the charged scalars 

i^hi i^k < niax(mfc,m/j) < lO^K^TeV. (60) 

On the other hand, if we use ^ 0.4(mfc/TeV), coming from r — > 3/i and |/^^| < 

0.1(m/i/TeV), coming from — > 67 we immediately obtain a lower bound on the masses of 
the scalars 

0.51 

mfc, iTLh > mm{mkmh) > TeV . (61) 



If we only use the r — > 3/i constraint in ( !58l) we find a bound on the |/^^| coupling 



/.|/,.r>2.6xl0~^( . /"^ (62) 
and using that ruh > min(mfc, mh) we find an absolute limit on the coupling 

\f^^r\>^^■ (63) 



Thus, using either the experimental bounds and/or the perturbativity bounds, we can also 
set upper and lower limits on the different couplings \g^^\ , /i and the interesting observables 
BR{fi — > 67) and BR{t —>■ 3/i). As discussed in section [XTTl the LHC will be able to find the 
doubly charged scalar of the model k^^ as long as it is lighter than about 1 TeV, thus it is 
interesting to know what are the constraints on the parameters of the model if rrik < 1 TeV. 
Following the procedure described above one can also put strong limits on the parameters of 
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General case 


TTijt < 1 TeV 


51 

' TeV < nih, ruk < lO^K^TeV 


TeV < ruh, mk < 24k^/^ TeV 


1 

— TeV < // < 10^ TeV 

K 


0.26 TeV < //< KTeV 


0.051 , „ , 0.01 ^ , , ^ 

r- < V < K, < \g^^\ < K 


0.051 , , , 0.1 ^ , , ^ 


BR{ii ei) > 1.0 X lO'^V'^^^ 


BR{^i ^ 67) > 8 X 10"^V«^^ 


BR{t 3/i) > 1.5 X 10"^^/k^2 


Si?(r ^ 3;u) > 2 X IO'^^/k'^ 



Table IV: Normal Hierarchy analytical constraints: we assume uJie — and central values for mea- 
sured oscillation parameters. For other experimental information we use 90% CL. The dependence 
on the perturbativity constraints is encoded in the parameter k and explicitly displayed. 

the model adding this additional constraint. We collect all the limits we obtain in table HVl 
It is important to remark the assumptions we use to obtain these bounds: we assume that 
because the small electron mass, as compared with the tau lepton and muon masses, — 0. 
We also take central values for the measured oscillation parameters. For other experimental 
information, bounds on branching ratios of rare processes, we use 90% CL limits. Finally the 
dependence on the perturbativity constraints is encoded in the parameter k and explicitly 
displayed. 

2. IH case 

The same kind of bounds can be obtained for the IH case with a few remarks. In the IH 
case ffj_r is not the largest coupling among the f's, since | V| ^ v^sisl V| with sis small. 
Thus perturbativity bounds should be applied to /e^. In addition the best experimental 
limit is also on /g^, |/e^| < 0.1(m/i/TeV). Then, it is convenient to write the main equations 
in terms of fefi instead of V- Finally in the IH hierarchy case there is the possibility of 
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General case 


m,t < 1 TeV 


^ TeV < mfc, m/j < 274k^ TeV 


^ TeV < mh.mk < Hk'^'^ TeV 


36 

TeV < /i < 274^5 TeV 


0.9 TeV < /i < KTeV 


0.095 , „ , 0.036 ^ , . 

^ < \ jeiA < K, ^ < \g^^\ < K 


0.1 , , , 0.36 


BR{^i ^ 67) > 2 X 10~^V«:^^ 


^ 67) > 1 X 10'^'^ /k^ 


BR{t ^ 3/x) > 2 X lO^^V'^^^ 


BR{t ^ 3/z) > 3 X 10"V«^^ 


0.0007/^3 < s?3 < 0.02 


0.018/k < < 0.02 



Table V: Inverse Hierarchy analytical constraints. 



cancellations for = vr which allow for a slightly smaller cj^^. We have in this case 

.;3l9..ll/.„P>T.3xlO-=|^H!5^(!n^. (64) 

Then we can repeat essentially the same arguments used for the NH, together with the 
upper limit on s^3, 5^3 < 0.02, to obtain lower and upper limits on the masses of the scalars, 
mh,mk, on the coupling |/e^|, which is related to l/erl and |/^t|, on the coupling \g^^\, 
related to and IqttI, and on the trilinear coupling fi. In addition, since in the IH case 
there is a strong dependence on 5^3 we can also set a lower bound on it. As in the NH case 
we also give the corresponding limits one would find under the assumption that the double 
charged scalar k^~^ is found at the LHC and, therefore, has a mass smaller than 1 TeV. We 
summarize all the limits in table |Vl 

B. Numerical analysis 

The information we obtained above is very useful; however, to obtain it we have made 
use of different approximations: 
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a) We assumed that the tUee, cUe^i, ^er can be neglected in front of the other couphngs. This 

approximation is reasonable because these w's are proportional to the electron mass, 
ujf.i = rriegeimi, which is much smaller than the other two lepton masses. However, it 
could happen that, for some reason, the Qei couplings are much larger than the others. 
It is therefore important to perform a complete analysis without this assumption. 

b) We took central values for the measured oscillation parameters. 

c) In the analytical limits we only used data from neutrino oscillations and bounds from 

r ^ 3/i and /i — > 67 (or lepton/hadron universality in the IH case) together with 
the perturbativity constraints. As discussed in section IIIBI there are many more 
experimental constraints that can affect the results and should be taken into account. 

It is clear that the only way to analyze the model without those approximations is by means 
of an exhaustive numerical exploration of the parameter space of the model. The basic tool 
to achieve this goal will be the use of Monte Carlo (MC) techniques; however, because of the 
large number of independent parameters and their diverse relevance, straightforward appli- 
cation of MC techniques is not sufficiently efficient and thus some additional considerations 
and refinements will be required. 

The crudest MC exploration of the available parameter space would involve random gen- 
eration of complete sets of 17 independent basic parameters, calculation of the corresponding 
predictions for the observables and finally an acceptation/rejection process in terms of the 
agreement between those predictions and the appropriate experimental constraints. Be- 
side the large number of parameters to be considered, the relations among them previously 
discussed render such a crude approach almost hopeless. 

Realizing that not all observables play an equal role, that is, some of them are much 
more informative or constraining than others, we can go one further step in the use of 
simple MC techniques: instead of the simplest MC outlined above, we can construct a 
MC process devised to automatically produce mass matrices in agreement with neutrino 
oscillation experiments. 

Knowing the masses and mixing angles, if we were to reconstruct the mass matrix Ai,^ 
using experimental input, the only missing ingredients would be the Dirac phase 6, the 
Majorana phase (see eqs. ( 138] - H0l) ) and the poorly known mixing ^13, for which we only 
know it is small and ignore its exact value or even if it is zero. Equation (|T3|l gives the mass 
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matrix Ai^ in terms of the new parameters - the Yukawa couphngs, the scalar masses and 
the trihnear couphng we can thus try to fix some parameters by matching the extracted 
Ai,^ from oscillation data and the calculated one. This procedure achieves two goals: it 
guarantees that neutrino oscillations are adequately produced and it reduces the freedom 
in parameter space entering numerical study by trading some of the couplings by measured 
neutrino oscillation parameters. For each set {Ou, On, O23, Aa, '^s, <P} we thus obtain nu- 
merical values for the entries in A^,^. We will then use eqs. fHSj) or eqs. fj46|) to fix /g^ and fer 
in terms of f^r and the generated mixing angles. Then come eqs. fH7|l : these three complex 
relations involve the six complex couplings gab, the trilinear coupling and the scalar masses. 
Together with m^, ruh and fi, knowing three independent QabS in eqs. fHTI) will be sufficient 
to fix the remaining ones; effectively this means that we will generate Qee, gefj, and Qer, and 
thus automatically fix gf^^, g^r and Qrr- Notice that this phase convention is compatible 
with the standard choice for the neutrino mass matrix, eqs. fl37H38p . 

To summarize, by generating five quantities - Aa, As, 9ij - according to experimental knowl- 
edge, two phases - 6 and (j) -, one real coupling /^t, two masses - nih and rrik the trilinear 
coupling fi and three complex gab, we are spanning the 12 moduli and 5 phases needed to de- 
scribe the model. That is, instead of the crude and utterly inefficient Monte Carlo procedure 
in terms of {mk,mh, fi, fab, gab}, we can use {sfj,AA,As,S,(j),ff,r,mh,mk,fi,gee,gefigeT} to 
explore the whole parameter space and guarantee the agreement with neutrino oscillations 
results prior to the use of the remaining experimental constraints, which constitute the next 
step, as they are then applied to accept/reject "candidate points". Notice that we have not 
specified the generation process of the different quantities involved: some discussion will be 
addressed below, the details of the numerical generation are summarized in table I VII 

Despite being operative and useful, this refined MC procedure is not the last word as one 
can do better. For this purpose we resort to the use of Markov Chain driven Monte Carlo 
(MCMC) processes of the Metropolis type. 

We have discussed the benefits of a refined simple MC procedure with respect to the 
crudest one: the next (and final) step to complete the numerical toolkit we use is the 
rather straightforward conversion of this refined MC into a Metropolis-like simulation which 
provides the results to be discussed. This is largely beneficial as (1) the efficiency of the 
MCMC process is sufficient to produce a reliable and smooth output for the different subcases 
under study, (2) the refined MC gives a helpful check of the consistency of the whole process. 
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Parameter 


Value 


Shape 


Parameter 


Value 


Shape 


As 


(7.9 ±0.3) X lO-^eV^ 


Flat 


Aa 


(2.5 ±0.25) X 10^3 eV^ 


Flat 


sin^ 6*12 


0.30 it 0.03 


Flat 


sin^ 6*23 


0.50 ± 0.08 


Flat 


sin^ 


[10-7; 2 X 10-2] 


Log flat 








S 


[0;27r[ 


Flat 


(t> 


[0; 2tt[ 


Flat 




[102;109] GeV 


Log Flat 


ruk 


[102;109] GcV 


Log Flat 


f jJ^T 


[10-7; k] 


Log flat 


A* 


[1; 10^0] GeV 


Log flat 


\9ee\ 


[10-7; 


Log flat 


arg(5ee) 


[0;27r[ 


Flat 




[10-7; 


Log flat 


arg(c^e/.) 


[0;27r[ 


Flat 


\9er\ 


[10-7; 


Log flat 


arg(5er) 


[0;27r[ 


Flat 



Table VL Numerical values. 



Let us now discuss the remaining details concerning the simulations; notice that, even if 
in the following we refer to the generation of parameters, which is appropriate for the MC 
process, the corresponding feature when dealing with MCMC processes is not generation 
but in fact how they enter the stepwise acceptance function, however, to avoid essentially 
duplicated discussions we will just mention what concerns the plain MC case. The main 
idea that drives our election of shapes and ranges of the different parameters is the need to 
perform an adequate exploration of the available parameter space, in particular one has to 
ensure that the regions which can yield interesting signals like the production of scalars at 
the LHC or branching ratios of exotic processes close to present experimental bounds are 
properly studied. In particular: 

• Neutrino oscillations results, i.e. the squared masses differences A^, A5 and the mixing 
parameters sin2^i2, sin2^23, are generated with fiat distributions within a ±1.64(7 
range around the quoted experimental value (this corresponds to 90% confidence level 

or probability region for a gaussian-distributcd uncertainty of the measurement). For 
sin^ ^13, however, we only have an upper bound: to span a reasonable range of values 
it is generated through a logarithmically fiat distribution from the upper bound down 
to very small values, cut off at 10-7. 

• The Dirac and Majorana phases, 5 and 4>, are generated according to flat distributions 
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spanning the whole available range [0; 27r[. 

• Concerning the independent Yukawa couplings fab and gab, moduli are generated 
through distributions logarithmically fiat to explore values that could potentially span 
several orders of magnitude. The applied upper bounds correspond to the different 
naturalness/perturbative cases under consideration. The arguments, as the Dirac and 
Majorana phases, are generated through flat distributions over the complete [0;27r[ 
range. 

• The masses of the new scalar fields and nih are generated with logarithmically flat 
distributions reaching up to 10^ TeV and bounded below at ~ 100 GeV to incorporate 
LEP-motivated constraints. In any case, the precise upper bound is irrelevant as far 
as it is beyond the analytic bounds presented in tables IIVI and |Vl 

• The trilinear coupling /i is also generated with a distribution fiat in its logarithm and 
limited by the perturbativity requirement. 

• We apply the remaining experimental constraints presented in section IIIBI in a sharp 
(straightforward acceptance or rejection) way: the only acceptable predictions are the 
ones within the quoted 90% CL ranges/bounds. 

• The simulation described in the previous points allows a very wide range of scalar 
masses. However, as discussed, the most interesting case is when m^. < 1 TeV and, 
therefore, the k'^'^ can be discovered at the LHC Thus, we have performed an inde- 
pendent simulation requiring < 1 TeV. 

• All the simulations are done for both the NH and IH cases and for two values of the 
perturbativity constraint k, = 1 and k = 5. 

The arbitrariness in the choice of priors and their impact in the final results is always a 
concern in this type of analyses. Because of this we have used several priors. In the case 
of the neutrino oscillation parameters we have repeated the analysis fixing the parameters 
at the central values, taking gaussian distributions around central values and using the fiat 
distributions we have finally presented here. The differences are marginal and we chose to 
present results for flat distributions because the results are slightly more conservative. For 
other parameters we also tried plain flat priors, but, specially for parameters that range 
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in several orders of magnitude, logarithmically flat distributions span more efficiently the 
parameter space. We checked that the distributions obtained for the observables considered, 
for which we found analytical lower and upper limits, do not depend too much on the choice. 

At this point the machinery used to perform the announced numerical studies has been 
completely presented, however some comments on the nature and interpretation of the out- 
put it produces are in order. 

The model under study naturally "lives" in a parametric space of high dimensionality. 
The standard statistical arsenal offers two different approaches to reduce this high dimen- 
sional information and produce tolerably low dimensional - usually one or two dimensional 
- output: the frequentist and the bayesian frameworks. Very schematically: 

• Frequentists assign confidence levels to the marginalized output through the best fit 
achievable with the remaining parametric freedom. 

• Bayesians assign probability densities to the marginalized output through the integra- 
tion over the remaining parametric freedom of the likelihood (of data for the given 
parameters) times the prior distribution/weight of parameters (this is just Bayes con- 
ditional probability inversion formula at work). 

Beside the long standing quarrel existing among practitioners of one or the other approach, 
both, with their information reduction schemes, unavoidably present some drawbacks to- 
gether with their statistical merits. As we want sensitivity to the parameter space available 
for the model to work, the procedure we have followed might look quite bayesian. Being 
aware of the dependence on prior election and the imprecise nature of the details behind 
many constraints^^, we do not intend at all to try and produce would-be highly orthodox 
statistical results neither interpret them as if they were so, and thus we have chosen the 
numerical details of the simulations - that is both experimental constraints and priors - as 
stated above without any further qualm. 

Usually the available information is just a 90% CL range and little additional knowledge on the distribution 
originating this range is given. Moreover the perturbativity constraints, as clearly seen in tables ITVl and Ivl 
are determinant and, like all theoretical constraints, no obvious confidence levels or statistical significance 
can be assigned to them. 



35 



V. RESULTS 



In this section we collect the main results of the paper. 

First we would like to study the impact of the different assumptions and experimental 
data in the analysis. 

To illustrate the impact of low energy constraints (/i 67, r — > 3/i, ...) we perform 
an independent simulation with only neutrino data and another simulation including all 
constraints (in the case of IH and k = 5) and represent the resulting distribution^^ of for 
the two simulations. We represent with a dashed line the results of the simulation only with 
neutrino data and with a solid line the simulation with all present experiments included. 

IH, K = 5, I' oscillations & full simulation 




ruk (GeV) 

Figure 6: Impact of low energy constraints (/i — > 67, r — > 3/i, ...): distribution when only 
neutrino data is included (dashed) as compared with the case in which all experiments are included 
(solid). Displayed data correspond to the IH case and k = 5. 

It is clear from the figure that only neutrino data allow (even prefer) relatively low masses 
of the order of 1 TeV or below. However, when low energy experimental data is included 
the lower limit on the rrik is pushed to larger values. We have to remark that the shape of 
the curves basically reflects the volume of the parameter space (from the other parameters) 

In the following we obtain the distributions as five million point samples from a MCMC exploration of 
the parameter space as described in section HVBI 



36 



and that the tails can be rather long, thus regions below rrik < 1 TeV, as we will see by 
performing simulations with rrik < 1 TeV, are not completely forbidden. 

Similarly, to illustrate the impact of the perturbativity constraint, we represent in fig. [7] 
the rrih distribution from a simulation with n = 1 and another one with k = 5 (both in IH 
case with all experimental information included). 



Figure 7: Impact of the perturbativity constraints: ruh distribution for k = I (dashed) and k = 5 
(solid). IH case with all experimental data included. 

We confirm with this figure the scaling of the bounds with the perturbativity assumptions, 
encoded in the parameter k, obtained analytically. Thus, for smaller values of k the allowed 
range of rrih is much smaller. For k, = 1, the preferred region of rrih is in the range ~ 
1 — 100 TeV (although with long tails) while for k = 5 this range is enlarged to 1 — 10000 TeV. 
Although it cannot be appreciated in this figure the lower bound on the mass is also sensitive 
to shown analytically (see table IVT). 

From figure [6] it is clear that all present data allow a wide range of k^^ masses however, 
the k~^~^, as discussed in IIIH can only be discovered at the LHC if < 1 TeV. Thus, this 
is the really interesting region of parameters to be studied. To study this region we perform 
an independent simulation implementing (we present results for the IH case and k = 5) all 
present constraints but assuming, in addition, that the k~^~^ has been discovered at the LHC 
and therefore has a mass rrik < 1 TeV. In fig. [8] we present the BR{fi — > 67) distribution in 
the two cases, general case and < 1 TeV. We see the dramatic impact in BR{fi —>■ 67) 



IH full simulation, fc = 5 & fc = 1 




-- K =5 
K = 1 




m,. (GeV) 
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of the discovery of the k~^^ at the LHC While present data allow branching ratios in the 
range 10'^^ - 10"", if the is discovered at the LHC then BR{fi ^ 67) ~ lO^^^ - lO^^^ 
and, therefore, will be probed at the MEG experiment. 

IH ft = 5, full simulation &: < ITeV 




BR(^t ^ 67) 



Figure 8: Impact of the discovery of the k'^^ at the LHC (m^ < 1 TeV): BR{fj, 67) distribution 
for the general case, IH and k = 5 with all present experimental results, (dashed) and requiring in 
addition that has been seen at the LHC (m^ < ITeV) k = 5 (solid). 

Until now we have presented results only for the IH case. In general, as also seen in 
our approximate analytical results, we expect roughly similar results in the NH and the IH 
case, except for a few parameters and/or observables. In particular we mentioned that in 
the IH case there is a preference for the Majorana phase around cf) = n because in that case 
there is a cancellation in the neutrino mass formulas. This is confirmed by the numerical 
calculation: in fig. [9] we represent the distribution of cf) for both, the NH case (dashed line) 
and the IH case (solid line). The data is taken from a simulation with k = 5, including all 
present experimental constraints and requiring that < 1 TeV. The distribution in the 
NH case is practically fiat, while in the IH case it is highly peaked at = tt. 

We also expect large differences in the NH and IH cases for the parameter sin^^is. In 
fig. [To] we represent the sin^ 6'i3 distribution for the two cases, NH and IH (k = 5, full data 
and m/c < 1 TeV). While it is constant in the NH case, in the IH case it is highly peaked at 
the maximum values allowed by present data and, in fact, there is an absolute lower bound 
on it, sin^ ^13 ^ 2 x 10^'^, which is not so far from the present upper limit, sin^ ^^13 < 2 x 10^^. 
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K = 5, rrifc < 1 TeV, NH k IH 




(j) (Majorana) 



Figure 9: Differences between tlie NH and IH cases: distribution of the Majorana phase (j): dashed 



in the NH case and sohd in the IH case. All present experimental data included and k = 5. 



K = 5. jrifc < 1 TeV. NH & IH 

I I I I , , r— f 



-- NH 
— IH 




-3 -2 

10 5x10" 10 

sin^ 6*13 

Figure 10: Differences between the NH and IH cases: distribution of sin^ ^13 as in fig. [9l 

Finally, to illustrate another interesting difference between the two cases, NH and IH, 
we have represented in fig. [11] the distribution of {m^)ee, the relevant matrix element in the 
neutrinoless double beta decay experiments. As before we assume k = 5 and < 1 TeV, 
however, this assumption has little influence on the result since, as shown in eqs. fHSt - 
mi), {mu)ee is a function of only the neutrino masses and the mixing angles. Thus, the 
shape of the curves and their position is just a consequence of the fact that the model 
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predicts a massless neutrino. In any case, from the figure it is clear that the model predicts 
(m^)ee ~ 0.001 - 0.005 eV in the NH and (m^)ee ~ 0.01 - 0.06 eV in the IH case. 

K = 5, mk < 1 TeV, NH & IH 



-- NH 
— IH 



10"^ 10"^ 10"' 

{m^)ee (eV) 

Figure 11: Predictions for (mi^)ee in the NH (dashed) and the IH (solid) cases. All present experi- 
mental data included and k = 5. 

Now, for the most interesting observables, those which give some interesting constraints 
or good perspectives in future tests, we present two-dimensional contour plots of the cor- 
responding distributions. The density of points has been calculated using 50 bins in a 
logarithmic scale for each axis. Then, 10 contour lines equally spaced, ranging from the 
maximum density to 1/1000 of it, have been represented. Thus, the last contour region, 
painted with a lighter color, represents the region with a small density of points but which 
still contains some points. For each pair of observables we present two plots. On the left we 
present the distribution when all present experimental constraints are imposed. The values 
in the interesting region for the LHC (m^ < 1 TeV) are very low but not zero. Thus, to 
better study this region we present in a second plot (right) the results of a simulation im- 
posing the additional constraint < 1 TeV. All results are given for the most conservative 
perturbativity assumption (k = 5). Scaling for more restrictive assumptions can be inferred 
from tables [IV] and El We discuss the relevant plots for both the NH and IH cases. 
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A. Normal Hierarchy 

Here we consider correlations among observables in the normal hierarchy case. 




Figure 12: NH: BR{^ 67) vs ruk; (left) general case (right) assuming that the k'^^ is seen at 
the LHC (mfc < ITeV). 

The most interesting observables of the model are m^, because, if small enough, it will 
be accessible at the LHC, and the BR{fi 67), which will be probed at a precision of the 
order of 10~^^ in the MEG experiment. Thus in fig. [12] we display the joint BR{fi —>■ e'~f)-mk 
distribution. We observe a clear correlation between these two observables. Present data 
(left) seems to prefer ~ lOTeV and BR{fi — * 67) above 10~^^, but a large region of 
values is not excluded rrik ~ 10^ — 10* GeV and BR{fi 67) ~ 10^^^ — 10^^^, however 
if the doubly charged scalar, k~^~^, is discovered at the LHC (m^ < 1 TeV) the situation 
changes dramatically and the simulation shows that the preferred values are in the upper 
range^'' ^ 600 GeV and BR{^ 67) ^ 10""^^ and that k~^~^ masses below 200 GeV and 
BR{n — > 67) below 10~^^ are very difficult to obtain in the model. 

Since BR{fi 67) depends more explicitly on rrih than on it is interesting to study the 
correlation between BR{fi 67) and rrih- In fig.[T3]we depict the allowed region in the plane 
BR{fi — > e'y)-mh, on the left for the general case and on the right for the case < 1 TeV. 
We see a strong correlation specially in the general case. From the figure on the left we see 
that in the general case rrih can be in a very wide range of values rrih ~ 10^ — 10* GeV but 
the preferred values are rrih ~ 40 TeV. On the other hand, if < 1 TeV the allowed range 

For specific numbers we take 3 contours in the plots. 
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of m/i is much smaller, ruh ~ 10^ — 10^ GeV, and is shifted to the lower edge. It still allows 
a large range of masses not accessible at the LHC 



NH full K=5 



NH full m^<l TeV k=5 
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Figure 13: NH: BR{fj, 67) vs rufi. 

T —>■ 3fi is mediated by k^^ exchange and governed by couplings which are largely fixed by 
neutrino mass data. Thus in fig. [H] we present results for the allowed values in the plane 
BR{t 3/i)-mjt. As expected there is a strong correlation in the two cases considered. 
In the general case the preferred values of BR{t —>■ 3fi) are in the 10~^^ range, although 
values as small as 10~^^ are allowed. In the < 1 TeV case the preferred values are 
BR{t —>■ 3fi) ^ 10~^ which is not so far from present limits, BR[t —>■ 3/i) < 3.2 x 10^^ , 
but values like 10"^'^ are not completely excluded. 
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Figure 14: NH: BR{t 3/i) vs m^. 
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B. Inverted Hierarchy 

In this section we consider correlations among observables in the inverted hierarchy case. 




Figure 15: IH: BR[ii — > 67) vs m^. 

In fig. [15] we represent BR{fi 67) vs m^. On the left for the general case and on the right 
with the additional assumption that the k^^ has been seen at the LHC [rrik < 1 TeV). The 
plots are similar to the plots obtained in the NH case, although slightly more restrictive. The 
allowed region in the general case is rrik ~ 10^ — 10^ GeV and BR{fi — > 67) ~ 10^^^ — 10^^^ 
and higher density values occur for ~ 50 TeV and BR{fi 67) ~ 10^^^. However, 
if nature chooses a k~^^ light enough to be produced at the LHC the model is much more 
constrained: it predicts that rrik is relatively large (masses below 400 GeV are only marginally 
allowed and the preferred masses are above 800 GeV). In addition BR{fi 67) > 10^^^ 
and the preferred range is above 2 x 10~^^. 

Figure [11] is also similar to fig. [13] but slightly more restrictive. In the general case we find 
rrih ~ 10^ — 10^ GeV and preferred values rrih ~ 40 TeV. For rrik < 1 TeV the allowed range 
of rrih is rrih ~ 500 GeV — 70 TeV which will make its detection at the LHC problematic. 
The constraints on BR{t — 3/i) are also stronger in the IH case, fig. [16], than in the NH 
case. The allowed regions are similar but more restrictive. Thus we find that in the general 
case the preferred values of BR{t 3/i) are in the 10~^^ range, although values as small as 
10~^^ are allowed. In the < 1 TeV case the preferred values are BR{t 3fi) ^ 3 x 10~^ 
(to be compared with present limits BR{t 3/i) < 3.2 x 10~^), but values like 10~^° are 
not completely excluded. 

Finally, as has been shown analytically, in the IH case there is a lower bound on sin^ 613. 
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Figure 16: IH: BR{fj, — > ej) vs nih- 
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Figure 17: IH: BR{t 3fj,) vs nik- 



Thus in fig. [TS] we represent sin^ versus BR{fi 67). We observe that there is a strong 
tendency towards relatively large values of sin^ ^13 even in the general case. The preferred 
values are in the sin^ ^13 ~ 0.01 region although values below 10^'^ do not seem completely 
excluded (the present upper limit is sin^ 6'i3 < 0.02). In addition smaller values of sin^ 6'i3 
require larger values of BR{^ 67). If we also require that the k'^'^ can be discovered at 
the LHC we find the preferred values of the model are constrained to a region sin^ ^13 ^ 0.01 
and BR{fi — > 67) > 10^^^. We also see that values of sin^6'i3 below 0.005 and BR{ji ej) 
below 10^^^ are very unlikely in this case. Notice that mixings as small as sin^^^i3 ~ 0.005 
will be tested in a near future 66j and that the MEG experiment, which will start this year, 
will probe BR{fi 67) at the level of 10~^^. 
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Figure 18: IH: sin^ 6*13 vs BR{^ ej). 

VI. CONCLUSIONS 



In this section we would like to highlight some of the most relevant conclusions drawn 
from the previous analysis. 

The explanation of observed neutrino mixings in the model sets very strong constraints 
on the structure of the couplings of the singly charged scalars to fermions, in particular, in 
the NH case the couplings must satisfy fer — ffirf^ — fefi- To fix the absolute value we need 
another observable, in the NH case we find that the strongest bound comes from fi —>■ e'j 
and tells us that \fei\ ^ 0.05(m/i/TeV) and If^rl ~ 0.1(m/i/TeV). In the case of IH the 
couplings must satisfy fer — —feii and |/er| ~ 5|/^r|- Then, from lepton-hadron universality 
we find \fei\ ^ 0.1(m/i/TeV), and |/^^| ^ 0.02(m;i/TeV), bounds which are similar to the 
bounds obtained from /i ej. 

The structure of the couplings of the doubly charged scalar is also very constrained by 
neutrino masses and mixings. In the case of NH they must satisfy, to a good degree of 
precision, that l^f^T-l ~ |fl'MAt|("^^/''Tir) and {grrl ~ \9fj.fj,\{^fi/''T^T)'^- In the case of IH this 
relation does not need to be satisfied exactly because the electron couplings gei can be 
relevant. However, we have seen that in a large region of the parameter space this relation 
is also required. Then, the best constraint comes from fi^fi^fi^ which tells us that 

Is'mmI ~ 0.4(mfc/TeV), l^f^^l < 0.024(mfe/TeV), l^f^^l < 0.0015(mfe/TeV). The g^i couplings 
are not constrained by neutrino data but are constrained by low-energy processes which are 
summarized in tables [T] and IIIII 

We find that the neutrinoless double beta decay parameter (my)ee is strongly constrained 
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in the model. We find 0.001 eV < \{m^)ee\ < 0.004 eV in the NH case and 0.01 eV < 
\{^v)ee\ < 0.06 eV in the IH case. This is just a consequence of the measured neutrino 
masses and mixings and the particular structure of neutrino masses of the model which 
predicts a massless neutrino. 

If this model is the right explanation for neutrino masses and if < 1 TeV the LHC 
will produce more than 10 events in the 4lep channel (see section IIIip There could be 
some dilution of the signal because the k can also decay into tau leptons or into two singly 
charged scalars but in the case of NH this can only be relevant for m^, > 1 TeV or if the 
singly charged scalar is light enough, rrik > 2mh. In the IH case the dilution of the Alep 
signal is a bit larger, still, most of the parameter space with < 1 TeV will give more 
than 10 events in the 4lep channel as long as the 2h channel is not open. 

If more than 10 events are produced at the LHC in the 4lep channel we find that BR{fi 
67) > 10~^3 in both the NH and the IH cases. These values are precisely the sensitivity 
expected in the MEG experiment at the Paul Scherrer Institute (PSI) which will start to 



run soon 231]. In fact one goal of the experiment is to obtain a significant result before the 
start of the LHC experiments. If this goal is achieved and nothing is seen we can reverse the 
argument and claim that it will be very difficult to find the charged scalars of this model at 
the LHC. We find that if BR{fi e-f) < IQ-^^ then, > 900 GeV and rrih > 600 GeV in 
the NH case and that both scalar masses will be above the TeV in the case of IH. 

It is also important to remark that the photonic vector form factor in muon-electron 
conversion in nuclei is enhanced with respect to the tensor form factor due to logarithmic 
corrections of loops with doubly charged scalars. Thus, if the current precision in muon- 
electron conversion experiments is increased in the next years, there will be additional tests 
on the model. 

If the doubly charged scalars are light enough to be produced at the LHC there are also 
interesting contributions to rare tau decays. For instance, in the IH case we find that most 
of the parameter space lies in the region BR{t — > 3/i) ^ 3 x 10~^, The NH case allows for 
slightly smaller branching ratios BR{t 3^) ^ 10^^. These results have to be compared 
with the the present limit, BR{t —>■ 3/i) < 3.2 x 10~^ or the ranges that might be explored 



We assume that, once efficiencies are taken into account, this corresponds to 2 reconstructed events 
5^,1521. 
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in SuperB factories ~ 10-^-10-^° |67|,l68|. 



The model gives a negative contribution to the = [g^ — 2)/2 of the muon. This means 
it cannot explain a positive deviation from the SM. Conversely, the precise measurements 
of set interesting constraints on the parameters of the model. 

In general it is much more difficult to satisfy all the constraints in the IH case than in 
the NH case. This can be seen in the MC acceptance rate which is much lower in the IH 
case than in the NH case. In fact, we have seen that satisfying all the neutrino mass data 
in the IH case requires certain cancellations in the neutrino mass formulas which imply that 
the neutrino Majorana phase (p cannot be zero, actually, most of the parameter space lies 
in the region e**^ ~ —1. 

Another interesting feature of the IH case is the emergence of a lower bound on the 
^13 mixing. We find that, even in the general case, most of the allowed parameter space 
requires sin^ 6'i3 > 0.001, although values below 10~^ do not seem completely excluded. In 
addition, smaller values of sin^6'i3 require larger values of BR{^ 67). If we also require 
that the k'^'^ is seen at the LHC through the 4/ep channel we find that values of sin^ 613 
below 0.005 are very unlikely (the present upper limit is sin^ ^13 < 0.02 and mixings as small 
as sin^^i3 ~ 0.005 will be tested in a near future |6(| ). 

In short, the requirement that the model is able to explain the observed pattern of 
neutrino masses and mixings places very strong limits on the parameters of the model. Thus, 
if the doubly charged scalar of the model is seen at the LHC through the Alep channel, the 
model predicts large contributions to several low energy processes, fi 67, r 3/i, hadron- 
lepton universality tests, which should be within reach of the next round of experiments. In 
particular, the MEG, /i ^ 67 experiment 23|] should provide results at the required level 
of precision before the start of the LHC experiments. If the k^~^ is discovered at the LHC 
and MEG does not see anything, the model will be in serious trouble. In all the other cases 
the model can fit all the data. Moreover, if the k~^^ is discovered at the LHC and MEG 
sees /i — > 67 in the allowed region, the model will be a serious candidate to explain neutrino 
masses. 
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